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ABSTRACT 

In  this  study  of  the  elastic  stresses  In  thick,  strip-wound 
magnet  colls  experimental  data  are  compared  with  theoretical  stresses 
based  on  three  physical  models.  It  Is  shown  that  the  floating  shell 
model  predicts  radial  dependence  of  the  circumferential  stresses  pro- 
duced  by  the  J  x  B  body  forces  better  than  the  monolithic  disc  or 
magnetic  pressure  models. 

The  fundamental  equations  relating  the  magnetic  body  forces 
to  the  stresses  In  the  conductors  are  derived  from  first  principles. 
The  instrumentation  for  measuring  stresses  within  magnet  windings  Is 
described.  A  survey  of  several  recent  magnet  stress  studies  Is 
Included. 
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NOMENCLATURE 

Inside  radius  of  cot  I 
Outside  radius  of  coll 
Magnetic  Induction 

Magnetic  induction  at  geometric  center  of  coll 

Half  length  of  coll 

Capacitance 

Electric  displacement 

Unit  vector 

Modulus  of  elasticity;  also  electromotive  force 
Magnetic  force 

Coll  field  factor  defined  tn  Equation  (7) 

Magnetic  field  Intensity 
Current  per  turn  of  conductor 
Current  density  In  conductor 
Constant  defined  In  Equation  (32) 

Inductance  of  coll 

Absolute  value  of  the  slope  of  the  BZ/BQ  versus  y  plot 

Number  density  of  current  carrying  charges 

Number  of  turns  In  coll 

Magnetic  pressure 

Electric  charge 

Electrical  resistance 
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r  Radial  coordinate:  also  radius  at  some  point  In  coll 

t  Time 

u  Radial  displacement 

v  Velocity 

V  Conductor  volume 

Vc  Capacitor  charge  voltage 

y  Value  of  ratio  BZ/0O  at  y  =  I 

z  Distance  along  axis  from  geometric  center  of  coll; 

also  axial  coordinate 

a  Ratio  of  outer  to  Inner  coll  radius,  82/^1 

B  Ratio  of  coll  half  length  to  coll  inner  radius,  b/a( 

Y  Ratio  of  Intermediate  radius  to  coll  Inner  radius,  r/aj 

6  Incremental  change 

A  Incremental  change 

e  Unit  strain 

0  Circumferential  coordinate 

A  Coil  space  factor 

u  Linear  permeability 

u0  Permeability  of  free  space 

v  Po I sson ’ s  rat  I o 

a  Unit  stress,  also  electrical  conductivity 

t  Diffusion  time  constant 

$>  Magnetic  flux 

ti)  Radian  frequency 
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Subscripts 

r  Radial  direction 

7  Axial  direction 

0  Circumferential  direction 

Superscripts 

*  Denotes  plane  strain  condition 

Denotes  peak  value 
"*■  Denotes  vector  quantity 
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CHAPTER  1 

INTRODUCTION 

I .  STATEMENT  OF  PROBLEM 

Air  core  electromagnets  are  used  extensively  to  generate  the 
high- Intensity  magnetic  fields  required  In  many  current  research  pro¬ 
grams.  The  field  attainable  with  a  given  magnet  Is  limited  by  the 

-V 

J  x  B  body  forces  exerted  on  the  windings.  This  body  force  (force  per 
volume),  gives  rise  to  stresses  In  the  current  conductors  proportional 
to  the  square  of  the  field  intensity  and  may  cause  mechanical  break¬ 
down.  A  magnet  coll  may  fall  by  bursting  or  by  localized  yielding  and 
subsequent  electrical  shorting. 

A  survey  of  the  literature  reveals  that,  although  the  subject 
of  magnet  stresses  has  been  treated  in  detail  Cl  through  93*,  quite 
diverse  results  are  obtained  by  use  of  various  models  which  have  been 
assumed.  The  Intent  of  this  study  was  to  establish  the  most  accurate 
of  several  models  for  the  calculation  of  stress  distributions  within 
the  windings  of  a  widely  used  type  of  magnet  construction. 

The  particular  magnet  construction  considered  herein  Is  the 
spirally-wound,  cylindrical  coll  with  Inside  radius  a^,  outside  radius 
a2,  and  length  2b  as  depicted  In  Figure  I.  The  coil  is  assumed  to  be 
uniformly  wound  with  strip  or  foil  conductor.  For  generality  the  coll 

^Numbers  in  brackets  refer  to  similarly  numbered  references 
in  the  bibliography. 
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Figure  I.  Cylindrical  magnet  coil  geometry  and  coordinate  system. 
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dimensions  are  expressed  In  terms  of  the  nondlmenslonal  shape  param¬ 
eters:  (I)  a,  the  ratio  of  sy  a|>  *2)  S,  the  ratio  of  b  to  aj, 
and  (3)  y.  the  ratio  of  an  Intermediate  radius  r  to  B|. 

II.  LITERATURE  REVIEW 

In  a  recent  work,  Hord  (I lH  presented  a  summary  of  several 
nrevlons  studies  C2  through  5]  of  the  stresses  In  various  type  magnet 
coils  and  developed  formulae  for  the  stresses  In  cylindrical  strip- 
wound  colls  by  use  of  three  different  analytical  models.  The  three 
models  Hord  considered  are  (I)  monolithic  disc:  this  model  assumes  an 
Isotropic,  homogeneous  cylindrical  disc  subjected  to  the  magnetic  body 
forces;  (2)  floating  shell:  this  concept  assumes  each  conductor  turn 
to  be  a  self-supporting  and  unrestrained  thin  wall  cylinder:  and  (3) 
men net l c  nressuro-  the  stress  Is  assumed  that  Incurred  by  a  thick 
wall  cvllnder  subjected  to  Internal  pressure  equal  to  the  magnetic 
energy  dens i tv. 

For  the  homogeneous  disc  model,  formulae  were  developed  for 
the  axial,  radial,  and  circumferential  stress  components.  The  radial 
and  circumferential  stresses  were  derived  for  both  the  plane  stress 
and  plane  strain  conditions.  The  circumferential  stress  component 
only  was  formulated  for  the  floating  shell  and  thick  wall  cylinder 
models.  The  formulae  are  dependent  on  the  magnetic  induction  $  and 
are  applicable  only  for  elastic  deformations  of  the  conductor  mate¬ 
rial.  Symmetrical  flux  distribution  and  uniform  current  density  were 
assumed.  Hord’s  analysis  shows  that  the  axial  component  of  stress  Is 
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on  the  order  of  20  per  Cent  of  the  maximum  circumferential  stress.  The 
radial  stress  component  was  considered  negligible.  On  the  basis  of  the 
mode  of  failure  of  one  magnet  constructed  of  a  stack  of  strip-wound, 
epoxy-bound  coils,  Hord  concluded' that  (I)  the  homogeneous  disc  model 
analyzed  under  the  assumption  of  plane  stress  is  applicable  for  tightly- 
wound,  epoxy  impregnated  coils,  and  (2)  the  floating  shell  model  best 
fits  loosely  wound  unbonded  coils. 

Daniels  [3]  also  dealt  with  stresses  in  a  strip-wound  coil. 

By  treating  the  windings  as  a  homogeneous  anisotropic  material,  the 
Insu lating  space  between  turns  was  accounted  for.  It  was  concluded 
that  each  turn  is  In  equilibrium  under  its  own  body  forces  and  that 
axial  compression  does  not  affect  the  magnitude  of  the  radial  and 
circumferential  stress  components.  Daniels  also  considered  the  sta¬ 
bility  of  coil  windings  against  deformation  and  concluded  that  certain 
of  the  outside  turns  are  intrinsically  unstable  due  to  negative  body 
forces  in  that  region  and  that  in  the  absence  of  elastic  restraining 
forces,  the  outer  turns  may  collapse  upon  adjacent  turns  resulting  in 
coil  failure  by  radial ly  unsymmetrlc  distortion  of  the  windings. 

More  recently,  Gersdorf,  et  al.  [73  treated  plastic  as  well 
as  elastic  deformations  In  magnet  coil  windings.  The  anisotropic 
nature  of  the  windings  and  the  conditions  of  stress  compatibility  were 
considered.  Results  similar  to  those  of  Hord  D3  were  obtained  for 
the  elastic  stresses.  However.,  it  was  pointed  out  that  a  positive 
radial  stress  is  not  possible  In  loosely  wound  coils  because  It  would 
separate  the  windings  at  a  very  low  value.  To  account  for  this,  the 
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windings  were  divided  Into  two  parts.  For  the  Inner  part  the  radial 
stress  was  assumed  to  be  zero  so  that  essentially  the  floating  shell 
condition  was  realized.  For  the  outer  part  of  the  windings,  the  solu¬ 
tion  for  the  circumferential  stress  was  very  close  to  Hord's  plane 
stress  solution. 

According  to  Furth  [63,  the  axial  magnetic  field  of  an  air 
core  solenoid  behaves  like  a  two-dimensional  gas  In  that  it  exerts  a 
radial  pressure  equal  to  the  magnetic  energy  density.  Furth  treated 
long  solid-bodied  coils  of  helical  or  single  turn  construction  in  terms 
of  a  thick  wall  cylinder  with  this  magnetic  pressure  applied  to  the 
Inside  wall  surface.  Equations  for  both  the  radial  and  circumferential 
stresses  were  formulated.  For  thick  coils  the  maximum  tensl le  stress 
was  indicated  to  be  equal  to  the  magnetic  pressure,  it  was  implied  that 
the  maximum  magnetic  field  H„,  attainable  with  such  a  magnet  is  equal  to 
(2om/uQ>1^^  where  om  Is  the  maximum  permissible  stress  for  the  magnet 
material . 

Ignatchenko  and  Karpenko  C83  also  have  studied  the  problem  of 
stresses  in  solld-bodled  coils.  Using  the  theory  of  maximum  shear 
stresses  I n  the  th I ck  wa 1 1  cy 1 1 nder  mode  I ,  they  cone  I uded  that  the 
maximum  achievable  field  is  (om/u0)^. 

Several  of  these  theories  and  models  were  employed  in  an 
attempt  to  determine  the  maximum  field  attainable  without  damage  to  an 
existing  large,  pulse  operated  magnet  used  to  generate  the  B-field  for 
an  experimental  magnetohydrodynamic  accelerator  [iOD.  The  magnet  is 
constructed  of  two  stacks  of  thin  col  I  discs  wound  from  flat  rectangular 
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magnet  wire.  Due  to  significant  differences  in  the  results  obtained, 
this  study  was  made  in  an  attempt  to  estab I ish' the  most  accurate  model 
for  such  a  coli.  Data  obtained  from  a  small  experimental  strip-wound 
coil  Instrumented  with  strain  gages  and  pulsed  from  a  capacitor  bank 
energy  storage  system  are  compared  with  stresses  calculated  by  use  of 
the  models.  The  International  System  of  Units  is  employed  throughout 
except  for  the  measured  and  calculated  stresses  presented  In  Chapter 
IV.  These  are  expressed  in  Engineering  System  Units. 
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CHAPTER  II 

BODY  FORCE  ANALYSIS 

I.  MAGNETIC  FIELD  DISTRIBUTIONS  AND  BODY  FORCES 

It  is  necessary  to  develop  analytical  expressions  for  the 
magnetic  body  forces  which  generate  the  stresses  In  magnet  windings. 
Consider  the  force  acting  on  a  charge  q  moving  with  velocity  v  through 
a  magnetic  field  H.  The  total  magnetic  force  acting  on  the  charge  Is 
given  by  the  Lorentz  force  C 1 1 71  as 

t  =  qv  x  5  ,  ( I ) 

•f  4 

In  which  B  ■  uH  Is  the  magnetic  Induction.  The  total  force  6F  on  a 
small  conducting  volume  6V  in  which  the  charges  move  essentially  with 

4  ■+ 

velocity  v  and  experience  the  induction  B  is  given  by 

i?  »  n  iV  q  (v  x  5)  (2) 


or  I n  the  limit  by 

dF  a  llm  .  * 

dV  =  6V  0  Cnq(v  x  5)]  *  J  xB  (3) 

where  n  Is  the  density  of  the  current-carrying  charges  and  It  Is  noted 
that  t  =  nqv.  Thus,  In  order  to  determine  the  magnetic  body  force  at 
any  point  within  the  current  conductors  of  a  magnet,  it  is  necessary 
only  to  know  the  local  values  of  5  and  In  this  study  as  for  most 
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Invest I qat Ions  of  this  problem  In  the  literature,  steady-state  current 
density  and  steady-state  magnetic  fields  are  assumed.  That  Is,  It  Is 
assumed  that  |T|  Is  constant  and  that  $  Is  equal  for  all  values  of  time 
to  the  steady-state  value.  The  validity  of  these  assumptions  for  a 
pulsed  coll  Is  discussed  In  the  latter  part  of  this  chapter. 

Brown  and  others  Cl 2]  have  computed  and  tabulated  numerical 
values  for  the  steady-state  axial  and  radial  components  (Bz  and  Br, 
respectively)  for  thick,  finite-length  solenoids  carrying  uniform 
current  density.  Values  of  Bz  and  Br  are  expressed  In  terms  of  Bq,  the 
flux  density  at  the  geometric  center  of  the  coll,  and  are  tabulated  as 
functions  of  a,  6,  and  y  for  Integer  values  of  o  and  B.  Typical  dis¬ 
tributions  of  Bz/B0  and  Br/BQ  along  with  a  typical  flux  linkage  are 
superimposed  on  a  coi  I  cross  section  In  Figure  2.  It  Is  apparent  that 
Br  induces  axial  forces  on  the  current  conductors  and  gives  rise  to 
axial  compressive  stresses  oz.  The  radial  forces  Induced  by  Bz  give 
rise  to  radial  and  circumferential  stresses  which  tend  to  burst  the 
coll.  Data  from  Brown's  tabulations  are  plotted  In  Figures  3,  4,  5, 
and  6  and  show  the  variations  of  Bz  and  Br  with  coil  geometry. 

In  a  more  recent  work.  Brown  and  Flax  Cl  3D  have  computed  and 
tabulated  numerical  values  of  Hz  and  Hr  for  semi-infinite  solenoids 
with  zero  inner  radius.  Each  field  component  Is  expressed  non- 
dimensional  ty  and  tabulated  as  a  function  of  nondimensiona!  field-point 

t  • 

coordinates.  Axially  symmetric,  uniform,  azimuthal,  steady-state  current 
density  Is  assumed.  By  superposition  of  appropriate  semi-infinite 
solenoids,  local  values  of  Hz  and  Hr  and,  subsequently,  of  Bz  and  Br 
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Figure  2.  Plot  of  typical  distributions  of  the  axial  and  radial 
components  of  magnetic  Induction  for  a  cylindrical  coll  carrying 
uniform  current  density. 
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Nondimensional  Radius,  y 

Figure  3.  Plot  of  the  axial  component  of  magnetic  Induction  for 
several  axial  positions  In  a  cylindrical  coll  carrying  uniform 
current  density. 
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0  1  2  3  4  5 


Nondimensional  Radius,  y 

Figure  4.  Plot  of  the  axial  component  of  magnetic  Induction  for 
several  values  of  o  and  6  for  a  cylindrical  coll  carrying 
uniform  current  density. 
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Figure  5.  Plot  of  the  radial  component  of  magnetic  Induction  for 
several  axial  positions  In  a  cylindrical  coll  carrying 
uniform  current  density. 
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Figure  6.  Plot  of  the  maximum  radial  component  of  magnetic 
induction  for  several  values  of  a  and  $  for  a  cylindrical 
coll  carrying  uniform  current  density. 


13 


AEDC-TR-69-264 


may  be  determined  everywhere  within  any  finite  cylindrical  coll  carry- 
I nq  uniform  current  density.  These  tabulations  facilitate  the  deter¬ 
mination  of  field  components  in  colls  for  which  a  and  B  are  not  Integers 
and,  thus,  have  been  utilized  In  this  study. 

If  the  coil  depicted  In  Figure  2,  page  9,  Is  wound  with  N  turns 
and  carries  a  uniform  current  of  I  amperes  per  turn,  the  current  den¬ 
sity  at  any  point  is 


NI 


e* 


(4) 


J  =  2ba,xlct-l)  Be 

where  the  space  factor  X  Is  the  ratio  of  the  volume  of  current  conductor 
to  the  total  vol 
be  expressed  as 


B<  = 


in  which  the  term 


of  the  coil. 

In  Appendix  A  it  is  shown 

NI 

(Ml 

To  +  (a2  +-B2)  ^  1 

2b (a- 1 ) 

l 5  1 

£nli  +  (1  +  b2)|/2  J 

i  written 

as 

ICT6  NI 

WCa-H  F(q>S) 

i)  -  (2ii| 

i  r a 

D  n  1  mv 

+  (a2  +  02)1/2  *| 

11  l  5  J 

*n  I 

L  1 

+  Cl  +  e2>'/2  J 

(6) 


(7) 


Is  cal  led. the  field  factor.  Combination  of  Equations  (3),  (4),  and  (6) 
yields  for  any  point  within  the  conductor; 
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-► 

dF 

dV 


106  Bo  T-.  + 

'aprrc^y  I  ee  *  <£V  ®r  +  Bz  ez> 


Thus,  the  radial  force  component  Is 


(8) 


d  Fr  I06  Bz  Bq 
dV  ”  a | A  F(a, B) 


and  the  axial  force  component  Is 


III  .  106  Bp  (|0) 

dV  a(A  F(a,B) 

As  suggested  by  Hord  DD  and  as  Is  apparent  from  Figure  3, 
page  10,  In  any  plane  perpendicular  to  the  coil  axis  between  z  =  -b  and 
z  =  b,  Bz  within  the  col  I  windings  (I  -  y  -  a)  may  be  expressed  as 


Bz  *  Bq  Cy  +  m  Cl-y)], 


(II) 


where  y  Is  the  value  of  Bz/B0  at  r  =  a(,  and  m  is  the  absolute  value  of 
the  slope  of  Bz/B0.  In  Appendix  B  It  is  shown  that  an  approximation 
for  the  distribution  of  Br  throughout  the  windings  is 


Br  =  Bof  ml  +  m2  s^prlj"  J  * 

Thus,  Equation  (9)  becomes 

d  Fr  I06  B02 

"W  =  a, A  F(a,e)  Cy  +  m<l-Y): 
and  Equation  (10)  becomes 


(12) 


(13) 
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d  Fz 
dV  =  " 


106  B02 
Fjx  F(a,8) 


(14) 


Equations  (13)  and  (14)  are  used  In  Chapter  III  in  the  development  of 
formulae  for  the  magnet  stresses. 


1 1 .  COMMENTS  ON  THE  STEADY-STATE  ASSUMPTIONS  ' 


In  formulating  expressions  for  the  stresses  I.n  magnet  colls, 
uniform  current  density  and  steady-state  magnetic  fields  usually  are 
assumed.  Thus,  the  validity  of  comparing  theoretical  stresses  with 
those  measured  In  a  rapidly  pulsed  magnet  coil  must  be  questioned. 
Consider  a  cylindrical  magnet  coil  spirally  wound  with  several  turns  of 
wide  strip  conductor  and  energized  by  a  half-sine  pulse..  The  stresses 
which  can  be  measured  by  use  of  electrical  strain  gages  are  those 
which  exist  on  the  surface  of  the  conductors  at  various  radial  posi¬ 
tions.  For  a  pulsed  coil  the  convenient  reference  point  for  data 
reduction  purposes  Is  the  peak  of  the  strain  pulse.  Thus,  the  points 
In  question  are:  (I)  does  the  magnetic  field  completely  diffuse  Into, 
the  conductor  and  (2)  does  the  current  assume  a  uniform  distribution 
during  the  pulse?  For  a  multi-turn  coll  the  two  problems  are  related 
in  some  complex  manner  and  to  the  knowledge  of  the  writer  an  exact 
treatment  has  not  been  reported  In  the  literature.  Analysis  of  a 
simplified  model  will  lend  some  Insight  to  the  problem. 

Consider  one  turn  of  the  cylindrical  coll  to  be  a  thin 
cylindrical  shell.  If  the  radius  of  the  cylinder  Is  much  greater  than 
the  wall  thickness,  the  problem  of  diffusion  of  the  magnetic  field 
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through  the  shell  can  be  treated  In  terms  of  diffusion  In  a  thin  flat 
plate.  Maxwell’s  equations  for  a  linear  Isotropic  conductor  are: 


3? 

aft 

7  X  H  =  J  + 


-► 

V 

-*■ 

V 


(15) 

(16) 

(17) 

(18) 


Other  Important  relationships  are: 

5  =  yH  (19) 

and,  In  the  absence  of  material  motion, 


J  =  at  ,  (20) 

where  y  and  o  are,  respectively,  the  permeability  and  conductivity  of 
the  conductor.  Substituting  Equation  (20)  into  Equation  (15)  yields 


(21) 


Also,  substituting  Equation  (19)  into  Equation  (16)  and  neglecting  the 
displacement  current  term  yields 


l 

v 


-► 

7  X 


S-J  . 


(22) 
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Now,  taking  the  curl  of  Equation  (22)  and  substituting  Into  Equation 
(21)  gives  , 

i—  7  x  (7  x  B)  =  -  !§.  .  (23) 

yo  3t 

But,  v  x  (7  x  B)  =  7(7  •  B)  B  and  7*5=0;  thus.  Equation  (23) 

becomes 


uo 


(24) 


which  Is  the  diffusion  equation. 

Equation  (24)  has  been  solved  Cl  lD  for  the  diffusion  of  a 
one-dlmenslonal  magnetic  field  Into  a  parallel  conducting  plate  as 
Indicated  In  Figure  7.  The  length  and  width  of  the  plate  were  assumed 
to  be  large  enough  compared  with  the  thickness  d  to  consider  the 
plate  as  Infinitely  large  In  the  x-  and  y-dl rectlons. 

For  an  applied  step  In  B  the  solution  for  B*  was  shown  DO 

to  be 


where . 


n  odd 


sin 


nirz  ■ 
-T3 


■°n+ 


°n  = 


_2_2 
n  if 


pod 


(25) 


(26) 


The  time  constant  for  the  fundamental  component  (n=l)  Is  defined  as 
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Figure  7.  Transient  magnetic  field  diffusion  In  a  thin  flat  plate. 
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T  = 


pad 


(27) 


I  n 

The  higher  harmonics  damp  out  much  faster:  therefore,  the  fundamental 
time  constant  Is  called  the  diffusion  time  constant.  It  can  be  shown 
Til]  that  after  a  time  t  >  3t,  the  field  Is  almost  completely  diffused 
into  the  plate.  For  copper  p  =  pOJ  and  a  =  5.9  x  10  at  ?93°K:  thus. 


t  =  7.52  d2  . 


(28) 


For  the  experimental  maqnet,  d  =  8.1  x  I0-4  meters  and  t  is  approxi¬ 
mately  5  x  I0“6  seconds. 

As  Indicated  above,  the  analogy  between  this  thin  plate  model 
and  a  multi-turn,  strip-wound  coil  is  to  let  the  cross  section  of  the 
nlate  in  the  x-y  plane  represent  the  cross  section  of  one  conductor 
turn  in  the  coil.  The  magnetic  field  then  corresponds  to  the  coll 
axial  field.  Obviously,  the  effects  of  the  self  axial  field  of  the 
turn  as  well  as  the  effects  of  radial  fields  generated  bv  the  other 
turn?  must  also  be  considered. 

Now  consider  the  problem  of  transient  current  density  dis¬ 
tribution  for  the  flat  plate  model.  From  Eouation  (16)  it  easily 
follows  upon  setting  p  =  pQ  and  neglecting  displacement  current  that 

4  =  _L  3B*  -  B  V*  I  nuz  \  ""  +A 

pQ  dz~  ~  ~  4  cos  —  U  .  (29) 

n  odd  1 


20 


AEDC-TR  -69-264 


Note  that  Jy  Is  produced  by  the  diffusing  magnetic  flux  density  and 
goes  to  zero  as  Bx  reaches  steady-state. 

The  diffusion  problem  for  the  real  coll  requires  a  more 
general  solution  of  the  diffusion  equation.  However,  the  numerical 
value  of  the  diffusion  time  constant  for  the  simplified  model  is  assumed 
Indicative  of  the  time  to  reach’  steady-state.  Thus,  for  variations  In 
coll  excitation  whose  characteristic  times  are  much  longer  than  t,  It 

4  4 

Is  assumed  that  the  steady-state  values  of  J  and  B  are  applicable. 
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CHAPTER  111 

STRESS  ANALYSIS 

In  theoretical  analyses  of  the  stress  distributions  within 
magnet  coils,  several  analytical  models  have  been  assumed  Ci — 9U.  Three 
such  models  are  DD:  (I)  an  isotropic,  homogeneous  cylindrical  disc 
subjected  to  the  magnetic  body  forces,  (2)  a  floating  shell  concept 
which  assumes  each  conductor  to  be  a  self-supporting  and  unrestrained 
thin  wall  cylinder,  and  (3)  a  thick  wall  cylinder  subjected  to 
Internal  pressure  equal  to  the  magnetic  energy  density.  The  deriva¬ 
tions  of  the  stress  formulae  developed  by  Hord  CG  $or  these  three 
models  are  outlined  in  this  chapter. 

I.  MONOLITHIC  DISC 


Consider  the  stresses  acting  on  an  incremental  volume  element 
of  a  thin  disc  taken  from  a  circular  cylinder  of  constant  wall  thick¬ 
ness  as  depicted  in  Figure  8.  Timoshenko  CI4U  shows  that  If  the 
cylinder  has  axi symmetric  loading  and  the  incremental  element  Is 
subjected  to  a  radial  body  force  dFr/dV  as  Indicated  In  Figure  9,  the 
equation  of  radial  equilibrium  is 


d  or 

“  °r " r  sr 


d  Fr 

r_w~  • 


(30) 


Plane  Stress 

For  the  case  of  plane  stress,  o_  =  0;  therefore  o0  and  o 

4  or 

may  be  written  C I 5U  in  terms  of  radial  displacement  u  as  follows: 
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dr 


Figure  8.  Equilibrium  of  an  Incremental  volume  element  of  a 
circular  cylinder  subjected  to  axlsymmetrlc  body  forces. 
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Figure  9.  Forces  acting  on  an  Incremental  element  of  a  circular 
cylinder  subjected  to  ax  I  symmetric  body  forces. 


24 


AEDC-TR -60-264 


0 

6 


E 


and. 


°r  * 


(31a) 


(31b) 


where  E  Is  the  modulus  of  elasticity  and  v  is  Poisson's  ratio.  By  com¬ 
bining  Equations  (13),  (30),  and  (31)  the  equilibrium  equation  becomes 


where. 


d^u  I  du  u 


=  KCy  +  m  (I-y)D, 


(32) 


106  BQ2  <|  -  v2) 

a | X  E  F(a,B) 


Equation  (32)  has  the  general  solution 


Cy  o 

u  =  C | r  +  — — •  +Kr^ 

The  constants  C|  and  C2  may  be  determined  by  combining  Equations  (31b) 
and  (33)  and  applying  the  boundary  conditions  or  =  0  at  r  =  aj  and 
r  =  a2;  thus, 

K  Q,  a, 

Cl  - - 24 -  »  (34a) 


and. 


C2  =  - 


K  02 


24 


(34b) 
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where  Q(  =  j^{Q2(  l-v)  +  8y(2+v)  +  m(7+5v)}  and, 


Q?  =  " 


l-v 


„2  7 

-4—  {8y(2tv)<  l-o)+m  C<7+5v)  -  8a(2+v)  +  02  (9+3v)J)  . 
o-l 


Now,  by  subs+t+u+fng  Equations  (33)  and  (34)  into  Equations  (31)  the 
circumferential  and  radial  stress  distributions  for  the  condition  of 
plane  stress  may  be  written  as 


~  24\  p(q  g)  Q|  ( l+v)+Q2( l-v)  -^2  +  3m( I+3v)y^-8( l+2v) (y+m)y j ( 35) 


I06  Br 


I 


(  I 


a-  = 


I06 


24A  F( 


R  ^ 

—  < 

0,8)  1 


0|(ltv)-02(l-v)  3m(3+v)Y2-8(2tv)(y+m) 

Y2 


Y  . 


(36) 


Plane  Strain 

For  the  case  of  plane  strain  the  axial  strain  ez  is  assumed 
to  be  zero.  For  this  case  oQ  and  or  in  terms  of  displacement  are 
written  as  Cl 6]: 


and, 


°6 


E 

(l+v)( l-2v) 


J(l-v) 


(37a) 


'The  stresses  are  In  newtons/meter2  when  B0 
webers/meter2.  To  convert  to  psi  multiply  B02  by  145  x 


is  expressed 

io-s. 


In 
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r  (l+v)(l-2v) 


(37b) 


By  substituting  Equations  (37)  Into  Equation  (30)  the  equilibrium 
equation  In  terms  of  displacement  is  obtained.  After  solving  for  u 
and  substituting  back  Into  Equations  (37)  the  circumferential  and  radial 
stresses  for  plane  strain  are  written  as 

.  .  '°6  Bo2 

°e  2WTv-\)¥Ta 


Tirr  -  «S  ^  ♦ 


8Y(y+m)(l+v)-3my2(l+2v)|  (38) 


2v) 


I06  B02 

°r  "  24X  Iv-iJPia 


TeT  [q.  .  Q.  <2-il 


+  8Y(y+tn)(2-v)-3my,t(3-2v)  ,  (39) 


-2v)J  , 


where 


Q*  B  -8y(2-v)[-^r  +  I  J  +  »(7- 2v)[-rI  ( 


I 


ma 


a2  -  , 


|  8a  (2-v)  +  3a2(2v-3)]  t 


and 


05  = 


a 


2  (a2- I ) (2v-l ) 


8y(2-v>( l-a)+m  [(7-2v)  -  8a(2-v) 


-  3a2(2v-3)D  >  . 


Axial  Compression 

Now  consider  the  case  of  axial  compression  in  the  monolithic 
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disc.  Axial  equilibrium  of  the  Incremental  element  shown  In  Figure  9, 
page  24,  requires 


d  oz  dFz 
dz~  =  W 


Combining  Equations  (14)  and  (40)  yields 


dc2 
d z~ 


IQ6  BpZ 
F(a,S) 


m  ] .  +  m2  sin 


(40) 


(41) 


After  integration  of  Equation  (41)  and  invoking  the  boundary  condition 
<xz  =  0  at  z  =  b,  the  axial  stress  distribution  within  any  plane  per¬ 
pendicular  to  the  axis  of  the  coil  Is  given  by 


I06  B02  (b2-z2)  r 
°z  *  "  2b  a,X  F(o,8)  [ml 


(42) 


II.  FLOATING  SHELL  METHOD 


The  so-called  floating  shell  concept  D3  yields  a  very  simple 
solution  for  the  circumferential  stress  component.  Consider  a  strip- 
wound  coll  in  which  each  turn  is  assumed  to  be  a  self-supporting  and 
unrestrained  thln-walled  cylinder  as  depicted  In  Figure  10.  The  radial 
force  dFf  on  the  element  of  circumference  subtended  by  the  angle  d6 
-acts  at  an  angle  0  to  the  horizontal,  and  consequently  Its  vertical 
component  Is  dFr  sin  0.  Thus  equilibrium  Is  satisfied  by 

2  o0  dr  dz  =  J  dFr  sin  0.  (43) 

v 
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y 


Ogdrdz  agdrdz 


Figure  10,  Equilibrium  of.  a  differential  length  of  a  thin  walled 
cylinder  subjected  to  axlsymmetrlc  radial  'body  forces. 


29 


AEDC-TR-69-264 


But  from  Equation  (13), 


tO6  B02 

dFr  =  a(X  F(a,8) 


Cy+m  ( l-v).1dV 


(44) 


and  In  cylindrical  coordinates  dV  =  r  dr  d0  dz.  Regarding  dr  and  dz 
as  Increments  rather  than  differentials.  Equation  (44)  becomes 


I06  Bq2 

a8  =  5a, X  F(a,P) 


/ 


Cy+m  ( l-y)]r  sin  e  d0. 


(45) 


Integration  between 'the  limits  of  0  and  ir  yields  for  the  circumferential 
stress 


I06  B02 

CT0  =  xT^b')"  tCy+m( 1-y)].  (46) 

III.  MAGNETIC  PRESSURE  CONCEPT 

Since  the  dimensions  of  |B|^/2u  are  those  of  energy  density, 
and  therefore  of  pressure,  the  magnetic  forces  may  be  Interpreted  as 
resulting  from  differences  in  magnetic  pressure.  If  the  magnetic 
pressure  concept  Is  applied  to  Incremental  elements  In  the  conductor 
volume,  It  may  be  shown  that  body  forces  result  which  are  Identical 
with  the  J  x  5  Interpretation. 

An  approximate  approach  [1,5,6]  treats  the  magnet  coll  as  a 
thick  walled  pressure  vessel  subjected  to  an  Internal  magnetic  pressure 
Pm  concentrated  at  the  Inside  wall,  where 
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Pm  = 


z/b=0 

Y=l 


(47) 


This  neglec+s  the  pressure  everted  at  the  outer  radius  of  the  coil  arid. 
It  Is  to  be  noted  for  the  coll  geometries  considered  here  that 


B 


z/b=0 
Y=  I 


(48) 


as  evident  from  Figure  4,  page  II.  In  the  literature  B  in  Equation 
(47)  Is  taken  as  Bq  which  is  only  true  for  an  infinite  solenoid,  or 
"scaling  factors"  are  Induced  to  account  for  the  larger  induction  at 
Y  c  I  for  finite  colls  Cl 7].  Pursuing  a  method  similar  to  that  used  in 
deriving  the  plane  stress  formulae,  the  circumferential  stress  for  a 
thick-walled  cylinder  subjected  to  Internal  pressure  Pm  is  found  to  be 


i2  P, 


°e  = 


m 


_  2  _  _  2 
a2  a| 


Rearranging  Equation  (49)  yields  the  form 


°0 


BZ  z/b=0 

v=i 

2p0(a?  -  1 ) 


(49) 


(50) 
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CHAPTER  IV 

EXPERIMENTAL  RESULTS 

Several  experiments  were  conducted  with  a  cylindrical  strip- 
wound  magnet  in  order  to  determine  the  variation  of  the  circumferential 
stress  with  radial  position  at  the  coil  midplane.  The  magnet  was 
energized  by  a  half-sine  current  pulse  from  a  30,000  microfarad  energy 
storage  capacitor  bank.  During  each  experiment,  simultaneous  measure¬ 
ments  were  made  of  BQ  and  the  circumferential  strain  eQ  at  one  of 
several  positions  within  the  coil  winding. 

I.  MAGNET  DESCRIPTION 

The  experimental  magnet  was  constructed  of  electrolytic 
tough  pitch  copper  strip  conductor  wound  with  alternate  layers  of  mylar 
film  insulation  onto  a  hollow  nylon  cylinder.  Prior  to  winding  the 
copper  conductor  was  vacuum  annealed  to  a  dead  soft  condition.  Yield 
strength  and  Young's  modulus  of  the  annealed  copper  were  experimentally 
determined  to  be  7000  psi  and  15  x  fO6  psi,  respectively.  Poisson's 
ratio  for  copper  is  0.355.  Physical  characteristics  of  the  coil  are 
given  in  Table  I.  During  winding  a  resistance-type  strain  gage  was 
bonded  to  the  conductor  on  approximately  every  fourth  turn.  Each  gage 
was  oriented  to  sense  the  circumferential  strain  component  at  the  coil 
midplane.  This  particular  coil  configuration  was  arrived  at  by  con¬ 
sidering  (I)  the  magnitude  of  strain  which  could  be  induced  with  the 
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TABLE  I 

PHYSICAL  CHARACTERISTICS  OF  THE  EXPERIMENTAL 
MAGNET  COIL 


Conductor  Thickness 

(meters  x  I0-2) 

.081 

Number  of  Turns 

N 

46.7 

1 nductance 

L  (hys  x  I0“6) 

101 

Resistance 

R  (ohms  x  10“^) 

6.8 

Coll  Half-Length 

b  (meters  x  10"^) 

1.91 

Col  1  Inside  Radius 

aj  (meters  x  I0"2) 

1.91 

Coi 1  Outside  Radius 

32  (meters  x  I0“2) 

6.10 

Radial  Shape  Parameter 

a  (a2/0|) 

3.2 

Axial  Shape  Parameter 

0  (b/aj) 

1.0 

Space  Factor 

0.90 
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available  energy  supply,  (2)  the  maximum  voltage  which  could  be  applied 
without  danger  of  breakdown  of  the  strain  gage  insulation,  and  (3)  the 
necessity  of  a  fairly  uniform  magnetic  field  over  the  axial  space 
occupied  by  the  strain  gages  to  minimi2e  errors  in  the  strain  measure¬ 
ments. 

Theoretical  values  for  the  axial  and  radial  components  of 
magnetic  induction  for  the  coil  as  determined  from  the  tabulations  of 
Brown  and  Flax  Tiji  are  plotted  In  Figures  il  and  12.  From  Figure  II 
the  values  of  y  a^d  !m|  in  Equation  II  are  found  to  bo  1.134  and 
0.634,  respectively.  Points  from  the  curves  In  Figure  12  are  used  In 
Appendix  B  for  calculations  of  mj  and  m2  in  Equation  (12). 

i I .  ENERGY  SUPPLY 

The  magnet  energv  supply  consisted  of  a  bank  of  950  oil- 
filled  125  uf  energy  storaqe  capacitors  with  a  charoino  system  and 
switching  circuit  to  deliver  the  erergy  to  the  magnet.  The  capacitors 
were  connected  in  a  series-parallel  arrangement  to  yield  an  equivalent 
capacitance  C  of  approximately  30,000  yf.  A  schematic  diagram  of  the 
system  is  shown  in  Figure  13.  The  switching  ignitron  was  a  National 
Electronics,  Incorporated  Model  NL-1037.  The  ignitron  was  fired  by  a 
i  kv  pulse  from  a  high  voltage  supply  one  millisecond  after  initiation 
of  the  data  recording  instruments.  During  the  experiments  the  capacitor 
charge  voltage  VQ  was  limited  to  500  volts  to  prevent  possible  break¬ 
down  of  the  strain  gage  insulation  and  subsequent  damage  to  measuring 
instruments.  Five  hundred  volts  was  not  sufficient  to  maintain 
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Figure  II.  Plot  of  the  axial  component  of  magnetic  Induction 
at  the  ml dp  lane  of  the  experimental  coll. 
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Figure  13.  Schematic  diagram  of  the  magnet  energy  supply  system. 
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Ionization  of  the  Ignltron;  therefore,  the  circuit  was  Interrupted  at 
zero  current  and  the  resulting  magnet  current  was  a  half-sine  pulse. 
Resistance  of  the  coll  is  negligible;  thus,  peak  value  of  the  current 
was  VQ  (C/L)*/2  (approximately  8000  amperes  for  VQ  =  500  volts)  and 
the  pulse  width  was  ir(LC)*^  (approximately  5.8  ml  1 1  (seconds) .  The 
high  voltage  relay  and  shorting  resistor  Indicated  In  Figure  13  were 
used  to  remove  residual  charge  from  the  capacitors  after  each  ’•shot". 

III.  INSTRUMENTATION 

Measurement  of  strain  In  the  magnet  conductors  was  accomplished 
with  a  carrier  amplifier  system  used  to  excite  and  measure  the  output 
of  a  strain  gage  bridge.  Magnetic  Induction  was  measured  by  Integrating 
the  voltage  output  of  a  search  coll  positioned  at  the  geometric  center 
of  the  magnet.  This  Integrated  signal  then  was  displayed  simultaneously 
with  the  strain  signal  on  a  dual  beam  oscilloscope.  A  camera  was  used 
to  obtain  permanent  records  of  the  oscilloscope  traces.  Additional 
equipment  Included  an  active  low  pass  filter  to  eliminate  high  frequency 
noise  generated  by  the  carrier  system  and  a  delay  network  for  synchro¬ 
nization  of  the  oscilloscope  trigger  with  firing  of  the  Ignltron.  A 
block  diagram  of  the  system  Is  shown  In  Figure  14.  A  more  thorough 
description  of  the  Instrumentation  Is  given  In  Appendix  C. 

IV.  DATA 

The  raw  data  obtained  consist  of  a  series  of  photographs  of 
oscilloscope  traces  of  the  strain  and  magnetic  induction  signals.  A 
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Figure  14.  Block  diagram  of  the  magnet  Instrumentation  system. 
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typical  oscilloscope  trace  of  the  signals  is  shown  In  Figure  15.  The 
top  trace  is  the  strain  signal.  By  use  of  the  appropriate  calibration 
factors  the  peak  values  of  the  oscillograph  deflections  are  converted 
to  circumferential  strain  e0  and  magnetic  induction  BQ.  The  stress 
then  Is  computed  by  use  of  the  relation  o0  =  e0E.  The  stress  may  be 

O 

expressed  In  a  rationalized  form  by  dividing  ofi  by  BQ  .  The  measured 
values  of  for  several  radial  positions  at  the  coil  midplane  along 

O. 

with  corresponding  values  of  <j0,  Bq,  and  c0/Bo  are  presented  in  Table 
II.  In  Figure  16,  experimental  values  of  o0/Bo^  are  plotted  versus  the 
nondlmensional  radius  y  to  show  the  variation  of  o0  across  the  radial 
dimension  of  the  magnet  winding.  Also  plotted  in  Figure  16  are  Equa¬ 
tions  (35),  (36),  (38),  (39),  (42),  (46),  and  (50'  as  evaluated  for  the 
test  magnet  configuration. 

V.  DISCUSSION 

As  can  be  seen  from  Figure  16  the  experimentally  determined 
stress  distribution  does  not  agree  closely  with  any  of  those  computed 
by  use  of  analytical  models.  However,  Equation  (46)  developed  from  the 
floating  shell  concept  predicts  a  distribution  which  has  the  same 
general  shape  as  the  experimental  distribution.  The  experimental 
values  of  a0  are  less  than  floating  shell  theory  near  the  Inside  radius 
and  are  maximum  at  a  greater  radial  position  than  the  calculated 
maximum.  Of  particular  Interest  is  the  fact  that  a  compressive  cir¬ 
cumferential  stress  as  predicted  by  Equation  (46)  was  measured  near  the 
outside  radius.  As  Indicated  by  the  curves  for  Equations  (35)  and  (38), 
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Figure  15.  Oscillograph  record  of  typical  strain  and  Induction 
signals. 
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TABLE  II 


TABULATED  VALUES  OF  e0,  o  ,  B  AND  o0/Bo" 
FOR  THE  EXPER I MENT AL0MAGNET  COIL 


Nondimensional 
Radius,  y 

£0 

(micro- inches/ inch) 

°0 

(psi ) 

(webers?meter2) 

o0/Bo2 

(ps 1 Cweber/meter2]2) 

1.04 

31.6 

474 

2.91 

56 

59.0 

885 

3.06 

95 

68. 1 

1020 

3.38 

89 

1.26 

63.6 

954 

3.05 

103 

90.8 

1360 

3.26 

128 

97.3 

1460 

3.26 

137 

84.2 

1260 

3.03 

137 

1.45 

76.0 

1  140 

2.87 

138 

84.2 

1260 

3.01 

139 

63.1 

947 

2.56 

144 

84.2 

1260 

2.85 

155 

1.62 

121.0 

1820 

3.26 

171 

105.0 

1580 

3.03 

172 

136.0 

2040 

3.44 

172 

104.0 

1560 

2.96 

178 

1.91 

121.0 

1820 

3.50 

149 

90.8 

1360 

2.96 

155 

132.0 

1980 

3.44 

167 
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TABLE  M  (continued) 


Nondlmensional 
Radius,  y 

60 

(ml cro- I nches/ I nch ) 

°0 

(psi ) 

( webers/me ter^ ) 

O0/Bo2 

( ps I Cweber/me+er2]2) 

2.06 

105.0 

1580 

2.97 

179 

1 1 1 .0 

1670 

2.97 

189 

150.0 

2250 

3.44 

190 

84.2 

1260 

2.56 

192 

2.29 

109.0 

1640 

3.08 

173 

124.0 

I860 

3.26 

175 

109.0 

1640 

3.03 

179 

1  13.0 

1700 

3.03 

185 

2.46 

60.5 

908 

3.03 

99 

34.2 

513 

2.21 

105 

63.6 

954 

2.91 

113 

68.1 

1020 

2.94 

1  18 

2.67 

45.4 

681 

3.44 

58 

42.1 

632 

3.03 

69 

45.4 

681 

2.98 

77 

26.3 

395 

2.21 

81 

2.83 

13.6 

204 

2.89 

25 

27.2 

408 

3.38 

36 

22.7 

341 

3.03 

37 

21.0 

315 

2.85 

39 
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TABLE  II  (continued) 


Nond i mens  1 ona 1 
Radius,  y 

Ee 

(micro- Inches/Inch) 

°e 

(psi) 

(webers/meter^) 

( p  s  i  JI  we  b  e  r /mete  r  2D  2 ) 

3.06 

-34.1 

-512 

3.03 

-56 

-36.3 

-545 

2.91 

-64 

-21.0 

-315 

2.10 

-71 

-43.1 

-647 

2.96 

-74 
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Rationalized  Stress,  c/B  -  psi/ (weber/meter 
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Nondiraensional  Radius,  y 

Figure  16.  Plot  of  the  theoretical  and  experimental  stress  distributions 
at  the  midplane  of  the  experimental  coil. 
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the  plane  stress  and  plane  strain  theories  predict  maximum  circumfer¬ 
ential  stresses  greater  than  either  the  measured  stress  or  floating 
shell  theory.  Stresses  calculated  by  use  of  the  magnetic  pressure 
method  (Equation  50)  are  lower  than  those  indicated  by  all  other 
methods.  Although  no  experimental  values  of  the  radial  and  axial 
stresses  were  obtained.  Equations  (36),  (39),  and  (42)  are  plotted  In 
Figure  16  to  Indicate  the  relative  magnitude  of  these  components  as 
compared  to  the  circumferential  component,  it  Is  clearly  seen  that 
the  circumferential  stress  is  predominant. 

Since  the  stresses  plotted  in  Figure  16  are  In  terms  of  BQ^, 
from  a  plot  such  as  this  for  a  particular  col  I,  the  nominal  maximum 
field  attainable  without  yielding  can  be  determined  for  a  given  con¬ 
ductor.  In  Table  ill  the  conductor  yield  strength  and  maximum  B0  are 

tabulated  for  several  conductor  materials.  The  tabulated  values  of 

2 

maximum  BQ  are  based  on  c/BQ  =  200.  Although  a  Bq  of  approximately 
2 

15  webers/meter  is  Indicated  for  the  hard  copper  and  copper  alloys, 

strip-wound  coils  generally  are  limited  to  approximately  10  webers/ 

2 

meter  because  of  Insulation  fatigue  problems. 
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TABLE  III 

MAXIMUM  Bq  ATTAINABLE  WITH  EXPERIMENTAL  COIL 
GEOMETRY  FOR  SEVERAL  CONDUCTOR  MATERIALS 


Conductor 

Material 

Nominal  Yield  Stress 
(psl ) 

Maximum  Bq 
(webers/meter2) 

Aluminum  (annealed) 

4,000 

4.5 

Copper  (annealed) 

8,000 

6.3 

Z 1 rcon 1 um-Copper 
(annealed) 

12,000 

7.7 

Beryl  Hum-Copper  No.  10 
(annealed) 

25,000 

11.2 

Copper  (Hard) 

40,000 

14.1 

Zirconium-Copper 

(hard) 

50,000 

15.8 
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CHAPTER  V 

CONCLUSIONS 

From  this  study  1+  must  be  concluded  that  of  the  analytical 
models  Investigated,  the  floating  shell  concept  most  nearly  predicts 
the  circumferential  stress  distribution  In  colls  loosely  wound  with 
strip  conductor.  The  plane  stress  theory  for  a  monolithic  disc  more 
closely  predicts  the  maximum  value  of  stress  but  not  at  the  actual 
position  of  occurrence.  The  monolithic  disc  model  probably  applies  for 
single  turn  or  helical  disc  type  coils.  The  magnetic  pressure  method, 
although  widely  used,  predicts  for  this  coll  a  maximum  stress  of 
approximately  one-half  that  measured  and  calculated  by  other  methods. 
The  measurement  of  a  negative  circumferential  stress  near  the  outside 
wall  of  the  coil  supports  the  conclusion  of  others  that  mechanical 
instability  can  be  anjroblem  In  strip-wound  colls.  The  difference 
between  the  floating  shell  theory  and  the  measured  stresses  may  be  a 
result  of  transient  effects  in  the  magnetic  field  and  current  density; 
however,  no  consistent  explanation  has  been  developed.  Another  pos¬ 
sible  source  of  error  is  the  effect  which  axial  forces  have  on  the 
circumferential  stresses.  These  effects  have  been  neglected.  The 
type  Instrumentation  employed  In  this  study  Is  used  routinely  and  Its 
performance  is  well  known.  -A  system  error  of  approximately  t  5  per 
cent  is  suspected.  The  greatest  error  Is  likely  to  occur  In  reduction 
of  data  from  the  oscilloscope  trace  photographs.  Other  possible 
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sources  of  Instrumentation  errors  are:  (I)  transient  response  of  the 
carrier  amplifier,  harmonic  filter  and  oscilloscope,  (2)  measurement  of 
the  shunt  calibration  resistor,  (3)  measurement  of  the  resistors  and 
capacitors  for  the  passive  Integrator,  and  (4)  signal  drift  during 
acquisition  of  the  calibration  traces.  To  the  knowledge  of  the  writer, 
measurement  of  stresses  within  electromagnet  windings  has  not  been 
attempted  before. 

Transient  diffusion  effects  certainly  require  further 
Investigation.  Of  special  Interest  should  be  the  transient  stresses 
at  the  conductor  surface  due  to  current  density  and  magnetic  field 
at  the  surface  upon  application  of  a  short  pulse  voltage. 
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APPENDIX  A 


DERIVATION  OF  THE  MAGNETIC  INDUCTION  AT  THE  GEOMETRIC 
CENTER  OF  A  THICK  CYLINDRICAL  MAGNET  COIL 


A  generalized  form  of  the  Blot-Savart  D8]  law  states  that 
the  magnetic  Induction  3  at  position  r2  due  to  a  current  I  at  position 
r(  Is  given  by 


B  <r2> 


da  x  ( r2— r | ) 


lr2-r| 


(A-l ) 


where  d£  Is  an  element  of  current  conductor  along  the  contour  C.  Con¬ 
sider  the  circular  current  filament  of  radius  a  as  depicted  In 
Figure  17.  If  the  magnetic  field  Is  to  be  calculated  at  a  point  on  the 
z-axls  the  following  relations  exist: 


da  =  ade  C-  sin  0  ey  +  cos  6  By),  (A-2) 

r2-r|  =  -  a  cos  0  ex  -  a  sin  0  e^  +  z  ez  . 

Substitution  of  Equations  <A-2)  Into  Equation  (A-l)  yields 

2ir 


3(z)  = 


*oI  f  < 

4’J  ' 


z  a  cos  0  ex  +  z  a  sin  e  ey  +  a2ez) 
(z2  +  z2)3/2 


de 


which,  after  integration,  becomes 
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Figure  17.  Axial  magnetic  field  of  a  circular  current  filament. 
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*,_v  _  yoI  a2  t 

B<z)  -  ez  • 


(A-3) 


Now,  apply  Equation  (A-3)  to  each  element  dz  of  a  cylindrical  current 
sheet  of  radius  a  and  length  2b  as  shown  In  Figure  18.  If  the 
sheet  Is  uniformly  wound  with  Nz  turns,  each  element  dz  contains 
Nzdz/2b  turns;  thus,  the  magnetic  Induction  at  point  zQ  Is  given  by 


B  (z  )  =  _° J^1  f  <|z  (A-4) 

z  o  2B  T  J  [(z0-z)2  +  a^p/2 

o 


dz 


Integration  of  Equation  (A-4)  yields 

yo  Nz1  /  sin.  02  “  sin  el 


BZ<z0)  - 


"ZF 


(A-5) 


where  0j  =  -  tan-1  (zQ/a)  and  02  =  tan-1  (2b  -  zQ)/a.  At  the  geometric 
center  of  the  coll  z0  =  b;  thus,  8|  =  -  tan-1  (b/a)  and  02  “  +an_l 
(b/a).  Making  these  substitutions  Into  Equation  (A-5)  yields  for  the 
axial  component  of  induction  at  the  geometric  center  of  the  current 
sheet 


Bz(0)  = 


Up  NZI  ^ 

2(a2+b2)l/2 


(A-6) 


Now,  consider  a  thick  coll  make  up  of  several  concentric 
current  sheets  as  depicted  In  Figure  19.  Each  element  da  contains 
Nrda/(a2-a|>  turns;  thus. 
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Figure  18.  Axial  magnetic  field  of  a  cylindrical  current  sheet. 


Figure  19.  Axial  magnetic  field  of  a  thick,  finite  length, 
cylindrical  coll. 
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or. 


a 


Bz(0) 


wo  NzNr  /  da 

2Ca2-a,  >  J  (a2+b.2)l/2  ' 

ai 


(A-7) 


Bz(0)  = 


Mp  NzNr 
2 ( a2-a | ) 


32  +  ^a2^+  b2) 1/2 
aj  +  C a 1 2+b^>  t /2 


(A-0) 


Now,  since  a2/a|  =  a,  b/aj  =  B,  yQ  =  4  x  I0“7,  and  N^p  =  N  =  total 
number  of  turns.  Equation  (A-7)  can  be  written  as 


B  (0)  =  £n 

z  a, (a- I ) 


a  +  (o2  t  B2) >^2 


I  t  (|  +  32) 1/2 


x  l(T7.  (A-9) 


The  magnetic  field  along  the  coll  axis  is  purely  axial;  therefore,  the 
magnetic  Induction  at  the  geometric  center  of  the  col  I  Is  written  as 


NI 


[2^Ll 

o+  <a2  +  S2)1'2 

m‘n 

*i  +  (i  +  e2)'/2 

x  I0"6.  (A-10) 
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APPENDIX  D 


DERIVATION  OF  AN  APPROXIMATE  EXPRESSION  FOR  THE  RADIAL 
COMPONENT  OF  MAGNETIC  INDUCTION  FOR  A  THICK 
CYLINDRICAL  MAGNETIC  COIL 


Since  the  radial  magnetic  field  of  a  thick  cylindrical  coll 
carrying  uniform  current  density  Is  approximately  sinusoidal  across 
the  radial  dimension,  Hord  CG  suggests  that  the  field  distribution 
can  be  approximated  by  an  expression  of  the  form 


Kq  +  (K  -  K0)sinir 
Mo 

where  for  each  point  along  the  coll  axis 


(B-l> 


Ko 


+ 


Y=o 


and  K  is  the  peak  value  of  Br/Bo.  Values  of  Br/B0  for  any  finite  coll 
have  been  computed  and  tabulated  Cl 2, 13].  By  plotting  several  values 
of  K0  and  Kq-K  versus  z/b  a  linear  approximation  can  be  made  by; 

Kq  =  m(  z/b,  R  -  Kq  =  m2  z/b,  (B-2) 

where  mj  Is  the  slope  of  the  Kq  versus  z/b  plot  and  m2  is  the  slope  of 
the  (K  -  Kq)  versus  z/b  plot  for  a  given  coil  geometry.  Substituting 
Equation  (B-2)  into  Equation  (B-l)  yields 

^■F  ["■♦  "2  »'»"  (£}]].  (B-3) 
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As  pointed  out  by  Hord  Cl3,  Equation  (B-3)  Is  valid  everywhere  within 
tne  volume  occupied  by  the  coll  winding;  also,  the  approximation  is 
closest  at  the  coll  end  planes  where  maximum  compressive  force  Is 
generated. 

* 

Figure  20  shows  a  plot  of  Kg  versus  z/b  and  K  -  Kg  versus 
z/b  for  the  experimental  coll  described  In  Chapter  IV.  Values  of 
KQ  and  K  -  KQ  were  determined  from  the  radial  field  plot  In  Figure 
12,  page  36.  From  Figure  20  values  of  m |  and  m2  are  found  to  be  0.285 
and  0.194,  respectively. 
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Nondimensional  Coil  Length,  z/h 


Figure  20.  Plot  of  Kq  and  (ft  Kq)  versus  z/b  for  the 
experimental  magnet. 
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APPENDIX  C 

INSTRUMENTATION 

In  order  to  determine  experimentally  the  distribution  of  the 
circumferential  stress  component  across  the  radial  dimension  of  the 
experimental  coll  several  resistance  type  strain  gages  were  bonded  to 
the  current  conductor  at  different  radial  positions  during  the  winding 
process.  The  gages  used  were  William  T.  Bean  type  EA-09-250BK- I OC 
which  Is  a  foil  gage  with  an  epoxy  film  backing.  The  bonding  agent 
was  Eastman  type  910  contact  cement.  Each  gage  was  mounted  astride  the 
coll  midplane  so  as  to  sense  the  maximum  circumferential  strain  at  that 
particular  radial  position.  A  small  cylindrical  search  coll  was 
positioned  coaxially  Inside  the  magnet  coll  and  astride  the  midplane 
to  sense  the  magnetic  Induction  at  the  geometric  center  of  the  magnet. 

A  Tektronix,  Incorporated  Model  502  dual  beam  oscilloscope 
equipped  with  a  Tektronix  Series  125  Land  type  camera  was  used  for 
recording  simultaneously  the  strain  and  the  magnetic  Induction  signals. 
Both  signals  were  half-sine  pulses  of  approximately  5.8  milliseconds 
duration. 


Strain  Measurements 

The  manner  In  which  the  resistance  of  a  strain  gage  varies 
with  applied  strain  Is  given  approximately  by  Poisson’s  ratio  as 


(G.F.), 


(C-l) 
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where  ARg  is  the  incremental  change  of  gage  resistance  Rg,  AL/L  is  the 
unit  strain  to  which  the  gage  is  subjected,  and  G,F.  is  the  gage 
factor  which  varies  with  the  type  of  gage.  Since  the  unit  strain 
e  =  AL/L,  Equation  <C-I)  can  be  written  as 


e  = 


ARg 


Rgifi.F.) 


<  C— 2) 


Thus,  by  measuring  the  change  in  resistance  of  the  gage,  the  strain  may 
be  determined. 

Consider  a  strain  gage  with  resistance  Rg  connected  in  one 
arm  of  a  Wheatstone  bridge  as  indicated  In  the  following  sketch. 
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If  the  bridge  Is  excited  with  a  voltage  Ej,  the  output  voltage  is 
given  by 


Ej  (RgR3  -  R2R4) 

0  =  (Rg+R2)(R3+R4)  * 


(C-3) 


Let  the  gage  be  subjected  to  a  strain  e  so  that  Rg 
Equation  (C-3)  then  becomes 


becomes  Rg  + 


Now, 


E 


o 


=  E 


If  Rg  =  R2  =  R3  =  R4  » 


(Rg+ARg)R3  -  R2R4 
( Rg+ARg+R2 i ( R3+R4 )  * 

.  a 

ARg,  then,  to  a  close  approximation. 


(C-4) 


Eo  =  T" 


Combining  Equations  (C-2)  and  (C-5)  yields 


(C-5) 


EQ  =  (e)  (G.F.)  .  (C-6) 

Thus,  the  output  voltage  of  the  bridge  is  directly  proportional  to  the 
strain  to  which  the  gage  is  subjected. 

In  order  to  discriminate  between  the  strain,  signal  and  the 
voltage  appearing  due  to  electromagnetic  Induction,  a  carrier  amplifier 
system  was  used  for  strain  measurement.  A  schematic  diagram  of  the 
system  is  shown  in  Figure  14,  page  39.  The  amplifier  used  was  a 
Consolidated  Electrodynamics  Corporation  (CEC)  Type  1-127  which  provides 
a  5  volt,  20  K  Hz,  carrier  signal  for  bridge  excitation.  In  the  CEC 
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Type  1-127  amplifier  the  output  of  the  demodulator  consists  of  the 
Information  signal  plus  the  harmonics  of  the  carrier  frequency  generated 
by  the  rectification  action  of  the  demodulator.  An  Internal  filter 
partially  suppresses  the  harmonics;  however,  it  was  necessary  to  use  an 
additional  filter  between  the  carrier  amplifier  and  the  oscilloscope 
In  order  to  obtain  a  desirable  signal  to  noise  ratio.  An  active, 
fourth  order,  Gaussian  to  6  db,  transitional  filter  having  a  2000  Hz 
half-power  point  was  used.  A  schematic  diagram  of  the  filter  Is  shown 
In  Figure  21.  A  Fourier  analysis  shows  that  all  components  of  a  five 
millisecond  half-sine  pulse  above  1000  Hz  are  negligible;  therefore, 
no  significant  loss  of  signal  can  be  attributed  to  the  filter.  • 


Magnetic  Induction  Measurements 


According  to  Faraday's  law  of  electromagnetic  induction  the 


electromotive  force  E  associated  with  a  change  In  magnetic  flux  $ 
through  a  circuit  of  N  turns  is  given  Cl8]  by 


By  definition 


(C-7) 


(C-8) 


where  S  Is  the  surface  bounded  by  the  circuit.  Therefore,  Equation 
(C-8)  can  be  written  as 
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E  = 


S»n  da 


(C-9) 


If  the  circuit  Is  rigid  and  stationary  the  time  derivative  can  be 
taken  Inside  the  Integral  where  It  becomes  a  partial  time  derivative; 
thus, 


E  o  -  N 


(C-IO) 


Since  this  must  be  true  for  all  surfaces  S,  It  follows  that 


E  -  -  NA  ,  (C-ll) 

where  A  Is  the  area  of  the  surface.  If  the  circuit  Is  located  In  a 
magnetic  field  such  that  Q  Is  a  function  of  time  only.  Equation  (C-ll) 
may  be  written  as 


E(t)  *  -  NA  .  (C- 1 2) 

at 

In  the  above  discussion  It  has  been  assumed  that  all  N  turns 
of  the  circuit  experience  an  equal  number  of  flux  linkages,  which  can¬ 
not  always  be  the  case.  However,  It  has  been  shown  C 1 93  that  for  a 
circuit  consisting  of  a  cylindrical  coll  of  length  2b,  inside  radius 
8|  and  outside  radius  +he  error  due  to  unequal  flux  linkages  Is 
negligible  If  the  following  relationship  exists: 
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y 


(C-13) 


Now  consider  the  RC  network  shown  In  the  following  sketch, 
R 


RL  eQ(t> 


Passive  Integrator 

The  frequency  domain  transfer  function  of  the  network  Is 


rL 

Eo  (g>)  =  R^l,  (C- 14) 

E"j  (u)  Jo)CRR|_ 

WTT[  +  1 


In  the  steady-state,  at  frequencies  for  which  jwCRRl/R  +  Rl  »  I,  Equa¬ 
tion  ( C— 1 4 )  may  be  expressed  as 


e0(+) 


(C-15) 


where 


If  the  Input  to  the  RC  network  Is  the  search  coll 
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output  voltage  given  by  Equation  <C- 1 2 > ,  the  network  output  voltage  Is 


B 


Integrating  and  solving  for  B<t)  yields 

$(t)  -  •  eQ<t>  .  (C-17) 

The  search  col  I  for  the  experimental  magnet  was  constructed  of  126 
turns  of  American  Wire  Gauge  Number  44  film  Insulated  copper  wire 
wound  on  a  cylindrical  bobbin.  A  balanced  differential  Integrator  was 
used  to  eliminate  common  mode  error  signals.  An  RC  time  constant  of 
0.1  second  was  selected.  The  output  of  the  Integrator  was  displayed 
on  the  oscilloscope  simultaneously  with  the  strain  signal.  The  Induc¬ 
tion  measuring  system  Is  shown  schematically  In  Figure  22. 

Data  Reduction 

Determination  of  numerical  values  of  strain  from  the  strain 
signal  photographs  was  accomplished  by  use  of  a  shunt  calibration  pro¬ 
cedure.  Consider  a  strain  gage  with  resistance  Rg  attached  to  a  member 
to  be  tested  and  connected  Into  a  measuring  system  as  Indicated  In  the 
following  sketch. 
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Search 

Coil 


Passive  Integrator  Oscilloscope 

Input 


Figure  22.  Schematic  diagram  of  system  for  measuring  magnetic 
Induction. 


70 


AEDC-TR-69-264 


Rc 

I — ww- 


-ww- 


^  Switch 


Measuri ng 
System 


Shunt  Calibration  System 

The  gage  Is  under  constant  or  zero  strain;  therefore,  the  resistance 
seen  by  the  circuit  Is  Rg.  If  a  calibrating  resistor  Rc  Is  shunted 
across  Rg,  the  equivalent  resistance  Re  presented  to  the  measuring 
circuit  Is 


=  *9^ 

Rg  +  Rc  * 

The  resistance  change  which  has  been  affected  is 


(C-18) 


1Rg  =  ^  %  Jhr  ■  <c-'s> 

g  c 

In  order  to  produce  such  a  resistance  change  In  terms  of  strain,  the 
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gage  mast  be  subjected  to  a  strain  e  such  that 

ARg  =  e(Rg)(G.F.)  (C-20> 

Therefore,  the  output  Indication  obtained  from  the  system  by  artificial 
means  corresponds  to  an  actual  strain  of 


*  ■  *  <Rgi(G.F.>  ’  <C-21 

Combining  Equations  (C-19)  and  (C-21)  results  In  an  expression  for  the 
simulated  strain  es  In  terms  of  known  parameters;  that  Is, 


= 


Rg  +  Rc  G*^* 


(C-22) 


By  double  exposure  of  the  film,  a  reproduction  of  the 
oscilloscope  traces  for  two  different  conditions  may  be  displayed  on 
one  photograph.  First,  with  the  bridge  excited  and  the  gage  subjected 
to  no  strain,  a  zero  trace  may  be  obtained.  Then,  by  shunting  the 
gage  with  a  calibration  resistor,  a  trace  corresponding  to  the 
simulated  strain  is  obtained.  The  displacement  ds  between  the  two 
traces  corresponds  to  the  simulated  strain  given  by  Equation  (C-22). 
Therefore,  the  actual  strain  from  a  data  photograph  ts  given  by 


x  k,  (C-23) 

where  d  Is  the  trace  deflection  due  to  actual  strain  and  k|  is  the 
oscilloscope  sensitivity  constant.  Corresponding  values  of  stress  were 
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computed  by  use  of  the  fundamental  relationship  between  stress  and 
strain  given  by 


a  =  e  E  ,  CC-24) 

where  E  is  Young's  modulus  of  elasticity  and  o  is  the  stress. 

Values  of  magnetic  induction  at  the  geometric  center  of  the 
magnet  coil  corresponding  to  peak  values  of  strain  were  determined  by 
use  of  the  relationship 


B0=^(d)(k2>,  (C-25) 

where  RC  Is  the  integrator  time  constant,  N  is  number  of  turns  in  the 
search  coil,  A  Is  the  effective  area  encompassed  by  the  search  coil,  d 
Is  the  signal  trace  deflection  on  the  photograph,  and  k2  is  the 
oscilloscope  sensitivity  constant. 
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